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r 7n 1 Abstract 

S^-T We study a supersymmetric version of the seesaw mechanism type-Ill. The model consists of 

the MSSM particle content plus three copies of 24 superfields. The fermionic part of the SU(2) 

_C triplet contained in the 24 is responsible for the type-III seesaw, which is used to explain the 

observed neutrino masses and mixings. Complete copies of 24 are introduced to maintain gauge 
coupling unification. These additional states change the beta functions of the gauge couplings above 
the seesaw scale. Using mSUGRA boundary conditions we calculate the resulting supersymmetric 
■ mass spectra at the electro- weak scale using full 2-loop renormalization group equations. We show 

that the resulting spectrum can be quite different compared to the usual mSUGRA spectrum. We 

T — l ' discuss how this might be used to obtain information on the seesaw scale from mass measurements. 

Constraints on the model space due to limits on lepton flavour violating decays are discussed. The 
main constraints come from the bounds on fi — > e^f but there are also regions where the decay 
r — > /i7 gives stronger constraints. We also calculate the regions allowed by the dark matter 
constraint. For the sake of completeness, we compare our results with those for the supersymmetric 
seesaw type-II and, to some extent, with type-I. 
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I. INTRODUCTION 



Supersymmetry offers a number of advantages compared to the standard model (SM). 
To name just a few, SUSY has a dark matter candidate, it can alleviate the gauge hierarchy 
problem and the minimal supersymmetric extension of the standard model (MSSM) leads 
to gauge coupling unification, if SUSY particles exist with masses of the order of the electro- 
weak scale. However, in the MSSM neutrino masses are zero, just as in the SM. Neutrino 
oscillation experiments , on the other hand, have shown that neutrinos have tiny, but 
non-zero, masses and that mixing in the leptonic sector is large 

From a theoretical point of view, if neutrinos are Majorana particles, all models of neu- 
trino mass at low energies reduce to the unique dimension-5 operator [6| 

{nfU= f -f{HL){HL). (1) 

Neutrino experiments determine only f a p/A, but contain no information about the origin of 
this operator, nor about the absolute size of A. If / is a coefficient 0(1), current neutrino 
data indicates A < O(10 15 ) GeV. This is the essence of the "seesaw" mechanism. 

One can show that there are exactly three different tree-level realizations of the seesaw 
mechanism Type-I is the well-known case of the exchange of a heavy fermionic singlet 
(sl-[lo|. Type-II corresponds to the exchange of a scalar SU(2) triplet fill . 12|. In seesaw 



type-Ill one adds (at least two) fermionic ££7(2) triplets to the field content of the SM jl3|. 
If in case of type-II and type-Ill models one would extend the usual MSSM by just the 
superfields responsible for neutrino masses and mixings, one would destroy the nice feature 
of gauge coupling unification as they belong to incomplete SU(5) representations. This 
problem is easily cured by embed ding the new states in complete £77(5) representations, 



e.g. in case of type-II in 15-plets [14j and in case of type III in 24-plets [15|. Note, that 
the 24-plet contains beside the £77(2) triplet also a singlet state which also contributes 
to neutrino physics and, thus, one has in this case actually a mixture between type-I and 
type-Ill. 

Understanding the nature of supersymmetry breaking by measuring the soft parameters 
will be one of the central tasks if signals of SUSY are found at the LHC. All the more 
so, since one can possibly gain some insight into the high energy scale physics from such 
measurements. Two kind of measurements containing indirect information about the seesaw 
scale in SUSY models exist in principle: lepton flavour violating (LFV) observables and 
sparticle masses. In case of seesaw t ype -I, low energy LFV decays such as li — > lj + 7 



an d lj — > 31 j have been calculated in |16M25|; \x — e conversion in nuclei has been studied 



26 



27] . The type-II model has received less attention, although it has actually fewer 
free parameters than type-I implying that ratios of LFV decays of leptons can act ually be 
predicted as a function of neutrino angles in mSUGRA, as has been shown in (3,0. A 
first study has been done in 29|. We stress that such a setup can not explain neutrino data 



unless non-renormalizable operators are added, as indeed is done in 29]. This is due to the 



need of generating a sufficiently large splitting between the Yukawa couplings of the singlet 
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and the triplets, which can not be obtained from RGE running only. Moreover, in the above 
publications for the type-II and type-Ill models only 1-loop RGEs have been used. However, 
we will show that using 2-loop RGEs is important for the calculation of the spectrum as 
this leads to a shift of the GUT scale. 

Measurements at colliders, once SUSY is discovered, can provide additional information. 
LFV decays of left sleptons within mSUGRA have been studied for type-I in 30|, |3l| and for 
type-II in 
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31 



Precise mass measurements, in particular of the sleptons and sneutrinos, 
might also show indirect effects of the seesaw 32|-|34j. As mentioned above, the additional 
heavy states of type-II and type-Ill lead to changes in the running of the beta functions 
and also of the mass parameters above the seesaw scale leading to changes of the spectrum 
at the electro-weak scale compared to the usual mSUGRA expectations. From different 
combinations of masses one can form "invariants", i.e. numbers which to leading order 
depend only on the seesaw scale 15|, although there are important corrections at 2-loop for 
the type-II 
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and, as we will show in this paper, also for type-Ill. It is also interesting to 
note, that the additional Yukawa couplings at the high scale can lead to a mass splitting 
between smuons and selectrons which in principle can be measured at the LHC: it has been 
shown in ref. (iB) that such a splitting may be constrained down to (9(10~ 4 ) for 30 fb^ 1 
of integrated luminosity. In mSUGRA, one expects this splitting to be tiny, whereas in 
mSUGRA plus seesaw significantly different masses are generated, as has been shown for 
type-I in ref. [36 . 



The modified spectrum also affects the calculation of the relic density. Assuming the 
standard thermal history of the early universe only four very specific regions in parameter 



space of mSUGRA can correctly explain the most recent WMAP data |37|. These are 
(i) the bulk region; (ii) the co-annihilation line; (iii) the "focus point" line and (iv) the 
"Higgs funnel" region. In the bulk, where the SUSY particles are relatively light, no specific 
relations among the sparticle masses exist. In the co- annihilation line the lightest scalar tau 
is nearly degenerate with the lightest neutralino, thus reducing the neutralino relic density 
with respect to naive expectations 38, Esj]. In the "focus point" line 39, 0] Q^oh 2 is small 
enough to explain QoMh 2 due to a rather small value of fi leading to an enhanced higgsino 
component in the lightest neutralino and thus an enhanced coupling to the Z-boson. Lastly, 
at large tan /3 an s-channel resonance pair annihilation of neutralinos through the CP-odd 



Higgs boson can become important. This is called the "Higgs funnel" region |41|. Also in the 
seesaw models of type-II and III these regions exist but the regions get shifted. Moreover, 
if the seesaw scale is sufficiently low the co-annihilation region disappears in type-II models 
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We will show that the same happens in case of the type-Ill model and we will contrast 
the results of this model with type-I and type-II models. 

The rest of this paper is organized as follows. In the next section, we first define the 
model. For completeness, and since we will compare the results for the different variants, we 
give the definitions for minimal type-I and type-II seesaws as well. We have used SARAH 



43H45J to calculate the full 2-loop RGEs, based on the general expressions given in [46| . 



We have, where possible, compared our results to previously available work and generally 
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found agreement. However, 47| have calculated 1-loop RGEs for all parameters and found 
some differences in case of the seesaw type-II to the RGEs published in 14j. Our calculation 
agrees with [47J. We then turn to the discussion of the resulting SUSY spectrum. The large 
changes in the spectrum affects the predictions for the rates of rare lepton decays, such as 
fi — > cy, and the relic density as discussed in section IHIi We present in section IIVI our 
conclusions. In the appendix we first summarize the procedure on how to obtain the RGEs 
for the soft SUSY breaking parameters from the beta functions and anomalous dimensions. 
We then give the formulas at the 1-loop and 2-loop level for these quantities for the seesaw 
models of type-II and type-Ill for an arbitrary number of new seesaw particles which are 
decomposed according to their SM gauge quantum numbers. 



II. MODELS AND SPECTRA 

In this section we briefly recall the main features of the three seesaw models. In models 
of type-II and III one adds particles charged under the SM gauge group. As they correspond 
to incomplete SU(5) representations, they would destroy the nice feature of gauge coupling 
unification. For this reason we add at the seesaw scale(s) additional particles to obtain 
complete SU(5) representations which we briefly review below. A more detailed discussion 
including the embedding in SU(5) models can be found in 4^]. 

In the subsequent sections we present the various superpotentials. In addition there will 
also be the corresponding soft SUSY terms which, however, reduce at the electro-weak scale 
to the MSSM one and, thus, are not discussed further. The additional terms of the soft SUSY 
breaking potential, due to the heavy particles, do not effect the discussion presented later on, 
as their effect is at most of the order Mewsb/ M seesaw and, thus, can be safely neglected. In 
this paper we will assume common soft SUSY breaking at the GUT-scale Mqut to specify 
the spectrum at the electro- weak scale: a common gaugino mass M1/2, a common scalar 
mass m and the trilinear coupling A which gets multiplied by the corresponding Yukawa 
couplings to obtain the trilinear couplings in the soft SUSY breaking Lagrangian. In addition 
the sign of the fi parameter is fixed, as is tan/3 = v u /vd at the electro-weak scale, where 
Vd and v u are the the vacuum expectation values (vevs) of the neutral component of Hd 
and H u , respectively. The models discussed below also contain new bilinear parameters in 
the superpotential leading to additional bilinear terms in the soft SUSY breaking potential 
which are proportional to B n of the MSSM Higgs sector. The corresponding RGEs decouple 
and their only effect is a small mass splitting between the new heavy scalar particles from 
the new heavy fermionic states of the order Bp I M. SP p . w ,,„ . This leads to a tiny effect in the 
calculation of the thresholds at the seesaw scale(s) (48| which, however, we can safely neglect. 
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A. Supersymmetric seesaw type-I 



In the case of seesaw type-I one postulates very heavy right-handed neutrinos yielding 
the following superpotential below Mqut'- 



Wi 

Wmssm 



w. 



MSSM 



U C Y U Q ■ H u - D c Y d Q ■ H d - E c Y e L ■ H d + fxH u ■ H d 
1 



N C Y V L ■ H„ 



-N C M R N C 



(2) 
(3) 

(4) 



where A ■ B = AiB 2 — A 2 B\ denotes the SU(2) invariant product of two SU(2) doublets. 
This model can be embedded in an SU(5) using the following SU(5) matter representations: 
1 = N c , 5m = {D C ,L} and 10m = {Q, U c , E c }. For the neutrino mass matrix one obtains 
the well-known formula 



rriv 



u YjM^Y u . 



(5) 



Being complex symmetric, th e lig ht Majorana neutrino mass matrix in eq. ([5]), is diagonalized 
by a unitary 3x3 matrix U 
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rh v = U T ■ m u ■ U . 
Inverting the seesaw equation, eq. (jSJ), allows to express Y v as (49( 



(6) 



(7) 



where the m v and Mr are diagonal matrices containing the corresponding eigenvalues. R is 
in general a complex orthogonal matrix. Note that, in the special case R = 1, Y u contains 
only "diagonal" products y/Miirii. For U we will use the standard form 



/ 



U 



C12C13 



S12C13 



si 3 e 



iS \ 



" S 12C23 — Ci2S23Si3e C12C23 — 512^235136 



S23 c 13 



X 



\ S12S23 — Ci2C23Si3e l<5 — C 12 S23 — ■5i2C 2 3S 13 e i<5 C23C13 / 



( e iai / 2 \ 
e iQ2 / 2 
1 



with Cij = cos 9ij and Sij = sin#y. The angles 812, 9 13 and ^23 are the solar neutrino angle, 
the reactor (or CHOOZ) angle and the atmospheric neutrino mixing angle, respectively. 5 
is the Dirac phase and ccj are Majorana phases. In the following we will set the latter to 
and consider in case of 5 only the cases and tt. 



B. Supersymmetric seesaw type-II 



In seesaw models of type-II one adds a scalar SU(2) triplet T to generate neutrino masses. 
As this triplet carries also hypercharge one has to embed it in a 15-plet of SU(5) which has 
under SU(3) x SU(2) x U(l) the following decomposition 14 



15 



S + T+Z 



(9) 



S ~ (6,1, 



2 
3 



T~ (1,3,1), 



(3,2. 



1, 
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One has to add two 15-plets 15 and 15 to avoid a chiral anomaly below the GUT-scale. The 
££7(5) invariant superpotential reads as 

W = -^Y 15 5 ■ 15 ■ 5 + -=Ai5 H ■ 15 • 5 H + ^A 2 5 H • 15 ■ 5 H + Y 5 10 • 5 • 5 H 
y/2 V2 V2 

+ Y 10 10 ■ 10 ■ 5 H + M 15 15 • 15 + M 5 5 H • 5^ (10) 

with 5h = {H c , H u ) and 5# = (H c , Hd). We do not show the part responsible for the ££7(5) 
breaking as we take the £77(5) only as a guideline to fix some of the boundary conditions 
at M GUT . Below Mqut in the ££7(5)-broken phase the superpotential reads 

W n = Wmssm + \{YtLT x L + Y s ]FSxI^ + Y z D c Z x L 

+ l={\ x HifA + X 2 H u f 2 H u ) + M T f x f 2 + M Z Z X Z 2 + M S SA (11) 

where fields with index 1 (2) originate from the 15-plet (15-plet). The second term in eq. ([111) 
is responsible for the generation of the neutrino masses yielding 

Note that 

Y T = U T ■ Y T ■ U , (13) 

i.e. Yt is diagonalized by the same matrix as m u . If all neutrino eigenvalues, angles and 
phases were known, Y? would be fixed up to an overall constant which can be easily estimated 
to be 

^ 10l5(3eV /oosevy (14) 



A 2 V m h 

In addition there are the couplings Y$ and Yz, which, in principle, are not determined by 
any low-energy data. In the calculation of LFV observables in supersymmetry both matrices, 
Yp and Y z , contribute. Having a GUT model in mind we require for the numerical discussion 
later the ££7(5) boundary conditions, apart from threshold corrections, Yt = Ys = Yz at 
Mqut- 

As long as Mz ~ Ms ~ Mt ~ M 15 gauge coupling unification will be maintained. The 
equality need not be exact for successful unification. In our numerical studies we have taken 
into account the different running of these mass parameters but we decouple them all at the 
scale Mt(Mt) because the differences are small. 

C. Supersymmetric seesaw type-Ill 

In the case of a seesaw model type-Ill one needs new fermions S at the high scale belonging 
to the adjoint representation of ££7(2). This has to be embedded in a 24-plet to obtain a 
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complete SU (5) representation. The superpotential of the unbroken SU (5) relevant for our 
discussion is 

W = V2l M Y^ M l H - ho M Y w 10 M 5 H + 5 H 2A M Y^ n 5 M + ^2A M M 24 2A M . (15) 

As above, we have not specified the Higgs sector responsible for the ££7(5) breaking. The 
new parts, which will give the seesaw mechanism, comes from the 24m- It decomposes under 
SU(3) x SU(2) x [7(1) as 

24 M = (1, 1, 0) + (8, 1, 0) + (1, 3, 0) + (3, 2, -5/6) + (3*, 2, 5/6) , (16) 
= B M + Gm + Wm + Xm + Xm ■ 

The fermionic components of (1,1,0) and (1, 3, 0) have exactly the same quantum numbers 
as v c and E. Thus, the 24^ always produces a combination of the type-I and type-Ill seesaw. 
In the £77(5) broken phase the superpotential becomes 

Wm = Wmssm + H U (W M Y N - J^B M Y B )L + H U X M Y X D C 

+-B M M B B M + \g m M g G m + ^W m M w Wm + X M M X X M (17) 

As before we use at the GUT scale the boundary condition Y^ = Y B = Yx and M B = Mq = 
My/ = Mx- Integrating out the heavy fields yields the following formula for the neutrino 
masses at the low scale: 

mv = "f {~T0 Y * M * 1Yb + \ Y w M wYw} ■ (18) 

As mentioned above there are two contributions stemming from the gauge singlet as well 
as from the £77(2) triplet. In this case the calculation of the Yukawa couplings in terms 
of a given high scale spectrum is more complicated than in the other two types of seesaw 
models. However, as we start from universal couplings and masses at Mqvt we find that 
at the seesaw scale one still has M B — My/ and Y B ~ Yw so that one can write in a good 
approximation 

m v = -v 2 u ^M^Y w (19) 

and one can use the corresponding decomposition for Yw as discussed in section III Al up to 
the overall factor 4/5. 



D. Effects of the heavy particles on the MSSM spectrum 



The appearance of charged particles at scales between the electro-weak scale and the 
GUT scale leads to changes in the beta functions of the gauge couplings (13 , [TH| . In the 
MSSM the corresponding values at 1-loop level are (bi, b 2 , 63) = (33/5, 1, —3). In case of one 
15-plet the additional contribution is Abi = 7/2 whereas in case of 24-plet it is Abi = 5. 
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FIG. 1: Mass parameters at Q = 1 TeV versus the seesaw scale for fixed high scale parameters 
mo = Mi/2 = 1 TeV, Aq = 0, tan /3 = 10 and fj, > 0. The full lines correspond to seesaw type-I, the 
dashed ones to type-II and the dash-dotted ones to type-Ill. In all cases a degenerate spectrum of 
the seesaw particles has been assumed. 



This results in case of type-II in a total shift of Abi = 7 for the minimal model and in 
case of type-Ill in Abi = 15 assuming 3 generations of 24-plets. This does not only change 
the evolution of the gauge couplings but also the evolution of the gaugino and scalar mass 
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15, ^8|. Additional effects on the 



parameters with profound implications on the spectrum 

spectrum of the scalars can be present if some of the Yukawa couplings get large 
which can also happen in type-I models (5l| . In Fig. [TJ we exemplify this by showing the 
values of selected mass parameters at Q — 1 TeV versus the seesaw scale for fixed high 
scale parameters m = M\/2 = 1 TeV and we have set the additional Yukawa couplings to 
zero. As expected, the effects in case of models of type-II and III are larger the smaller the 
corresponding seesaw-scale is. The scalar mass parameters shown are of the first generation 
and, thus, the results are nearly independent of tan/3 and Aq. For illustration we show in 
Fig. [2] the corresponding spectrum where we have fixed tan/3 = 10 and Aq = 0. 

We note that in all three model types the ratio of the gaugino mass parameters is nearly 
the same as in the usual mSUGRA scenarios but the ratios of the sfermion mass parameters 
change 15J, |28|. One can form four 'invariants' for which at least at the 1-loop level the de- 
pendence on Mi/2 and mo is rather weak, e.g. (rn^—m 2 E )/M^, (mQ—m 2 E )/Mf, (m 2 D —m 2 L )/Ml 
and (rriQ — mfy/Mf. Here one could replace M\ by any of the other two gaugino masses 
which simply would amount in an overall rescaling. In Fig. [3] we show these 'invariants' in 
the leading-log approximation at 1-loop order to demonstrate the principal behaviour for 
seesaw type-II with a pair of 15-plets and seesaw type-Ill with three 24-plets. From this one 
concludes that in principle one has a handle to obtain information on the seesaw scale for 
given assumptions on the underlying neutrino mass model, if universal boundary conditions 
are assumed. For the type-I, i.e. singlets only, of course Abi = and no change with respect 



to mSUGRA are expected. If, for example, the seesaw III model would be realized in nature 
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FIG. 2: Example of spectra at Q = 1 TeV versus the seesaw scale for fixed high scale parameters 
mo = Mi/2 = 1 TeV, tan/3 = 10 and a > 0. On left panel Mf l ,m^o.m-+ while on the right panel 
we have Ma, m v o, m~+. The line codes are as in Fig. [TJ 




Mig = M 24 [GeV] 



FIG. 3: Four different "invariant" combinations of soft masses versus the mass of the 15-plet or 24- 
plet, M15 = M24. The plot assumes that the Yukawa couplings are negligibly small. The calculation 
is at 1-loop order in the leading-log approximation. The lines running faster up towards smaller M 
are for type-Ill seesaw, the values for type-II seesaw are shown for comparison. 

with three 24-plets having similar masses around 10 13 GeV one could e.g. show that the 
corresponding ratios cannot be obtained with one pair of 15-plets in the seesaw II model, 
thus excluding this possibility. However, taking the seesaw II with two pairs of 15-plets one 
would obtain similar ratios as in this case the corresponding additional beta-functions at 
1-loop would be Abi = 14, e.g. nearly equal to our seesaw III model. 

The leading-log approximation gives only the general trend, but there is an important 
dependence on the SUSY point chosen. In Fig. H]we show as illustration (m 2 L — m E )/Mf 
and {ttlq — m 2 E )/M 2 for different mSUGRA points and at different loop orders: the dashed 
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M 24 [GeV] M 24 [GeV] 

FIG. 4: The limits of the invariants in seesaw type-Ill models. Left: {m\ — m^/Mf, right 
(mQ—rn^)/Mf. The blue lines are for SPS3, the light blue one for mo = 500 GeV and Af 1 / 2 = 1 TeV, 
and the red one for uiq = M^m = 1 TeV; full (dashed) lines are 2-loop (1-loop) results. The black 
line is the analytical approximation, for comparison. 



lines are at 1-loop level whereas the solid ones are at 2-loop level. The points considered are 
SPS3 with m = 90 GeV, M 1/2 = 400 GeV, A = 0, tan /3 = 10, [i > and for the same 



values of Aq and tan (3 two points with M%/2 = 1 TeV: ttlq = 500 GeV and mo = 1 TeV. The 
black line shows for comparison the leading-log approximation. We observe that usually the 
approximation gets worse for lower values of M24 and this is even stronger at the 2-loop level 
which is a consequence of the large coefficient in the beta functions at the 2-loop level, see 
e.g. appendix IA4I Nevertheless, one sees that in general it gives the correct trend, but it 
might even fail completely, e.g. in the case of Mi/ 2 = tuq — 1 TeV. The reason for the drop 
around M24 ~ 3.5 x 10 13 is that the difference between the parameters goes to zero as can 
also be seen from the right of Fig. [TJ see also discussion below. 

Last but not least we note that the use of the 2-loop RGEs leads to a shift of Mqut from 
about 2xl0 16 GeV for 24-plet mass of 10 16 GeV to about 4xl0 16 GeV for 24-plet mass of 
10 13 GeV, which is part of the differences between the 1-loop and 2-loop results in Fig. HI 
Here Mqut is defined as the scale where the electro-weak couplings meet, e.g. gu(i) — 9su(2)- 
This implies also that there is some difference for the strong coupling which is, however, 
in the order of 5-10% which can easily be accounted for by threshold effects of the new 
GUT particles, e.g. the missing members of the gauge fields and the Higgs fields responsible 



for the breaking of the GUT group [53j. A second reason why the deviations between the 
leading log calculation, the case of 1-loop and 2-loop RGEs gets larger for smaller seesaw 
scale is that the increase of the beta coefficients implies larger values of the gauge couplings 
at the GUT scale. This implies that one reaches a Landau pole for sufficiently low values 
of the seesaw scale. As an example we show in Fig. [5] the value of the gauge coupling at 
Mqut = 2 x 10 16 GeV as a function of the seesaw scale for type-II with a pair of 15-plets 
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FIG. 5: Values of the gauge coupling at Mqvt = 2 x 10 GeV as a function of the seesaw scale, 
black lines seesaw type-II and green lines seesaw type-Ill with three 24-plets with degenerate mass 
spectrum; full (dashed) lines are 2-loop (1-loop) results. For the calculation of the electroweak 



(black lines) and type-Ill with three degenerate 24-plets (green lines). In both cases the 
2-loop RGEs imply a larger gauge coupling for a fixed seesaw scale. One sees that in case of 
type-II (type-Ill) in principle one could reach a seesaw scale of about 10 8 GeV (10 13 GeV). 
However, we believe that we can no longer trust even the 2-loop calculation for such large 
values of the as the neglected higher order terms become more and more important. 
Especially, we should not trust the "turn-over" of the invariants in Fig. H] for very low values 
of the seesaw scale, since the numerical calculation at these points is already very close to 
breaking down. 

We would also like to mention that, in the numerical calculation we find very often that 
one of the scalar masses squared, in particular staus and/or sbottoms, gets large negative 
values already for values of the seesaw scale larger than the Landau pole and thus we can 
not go to values of the seesaw scale as low as the examples shown in Fig. H]in many SUSY 
points. 

E. Lepton flavour violation in the slepton sector 

>From a one-step integration of the RGEs one gets assuming mSUGRA boundary con- 
ditions a first rough estimate for the lepton flavour violating entries in the slepton mass 
parameters: 



threshold the spectrum corresponds to m,Q = M\/2 = 1 TeV, Aq = 0, tan/3 = 10 and fi > 0. 




(20) 



(21) 
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for % 7^ j in the basis where Y e is diagonal, = Itl(Mgut / 'Mi)5ij and Yft is the additional 
Yukawa coupling of the type-k seesaw at Mqvt (k = I, II, III). We obtain 

9 

aj = 1 , an = 6 and am = — . (22) 



Note, that in case of the type-II the matrix L is degenerate and thus can be factored out. All 
models have in common that they predict negligible flavour violation for the right-sleptons 

m 2 E>ij * 0. (23) 

We know that these approximations work well only in case of the type-I models. Nevertheless 
they give a rough idea on the relative size one has to expect for the rare lepton decays U — > Ijj 
which very roughly scale like 

Br{h -»■ Ijj) oc a 3 mf } m ^ tan 2 0. (24) 

where fh is the average of the SUSY masses involved in the loops. Note, that for a given 
set of high scale parameters both, the different size of the flavour mixing entries and the 
changed mass spectrum, play a role. 



III. NUMERICAL RESULTS 

In this section we present our numerical calculations. All results presented below have 
jeen obtained with the lepton flavour violating version of the program package SPheno 
The RGEs of the seesaw II and seesaw III models have been calculated with SARAH 



54,55 

43l445| . All seesaw parameters are defined at Mqjjt and as mentioned in the previous section 
we require for models of type-II the boundary condition Yz = Y$ = Yt and Mz = Ms = Mt 
and in case of type-Ill models Yj\r = Yb = Y\y and Mb = Mq = Mw = Mx- We evolve 
the RGEs to the scale(s) corresponding to the GUT scale values of the masses of the heavy 
particles. The RGE evolution implies also a splitting of the heavy masses. We therefore 
add at the corresponding scale the threshold effects due to the heavy particles to account 
for the different masses. In case of type-Ill models off-diagonal elements are induced in the 
mass matrices. This implies that one has to go the corresponding mass eigenbasis before 
calculating the threshold effects. We use 2-loop RGEs everywhere except stated otherwise. 
In the appendix we give the necessary ingredients on how to obtain them in the seesaw 
type-II and III models. The analogous anomalous dimensions for the type-I model can be 



found in 56 



Unless mentioned otherwise, we fit neutrino mass squared differences to their best fit 



values [5] and the angles to tri-bi- maximal (TBM) values |57l |. Our numerical procedure is 
as follows. Inverting the seesaw equation, see eqs. f fl2l) and (118jl . one can get a first guess 
of the Yukawa couplings for any fixed values of the light neutrino masses (and angles) as a 
function of the corresponding triplet mass for any fixed value of the couplings. This first 



12 



guess will not give the correct Yukawa couplings, since the neutrino masses and mixing 
angles are measured at low energy, whereas for the calculation of m u we need to insert the 
parameters at the high energy scale. However, we can use this first guess to run numerically 
the RGEs to obtain the exact neutrino masses and angles (at low energies) for these input 
parameters. The difference between the results obtained numerically and the input numbers 
can then be minimized in a simple iterative procedure until convergence is achieved. As 
long as neutrino Yukawas are VY^ < 1 we reach convergence in a few steps. However, in 
seesaw type-II and type-Ill the Yukawas run stronger than in seesaw type-I, so our initial 
guess can deviate sizable from the correct Yukawas, implying in general also more iterations 
until full convergence is reached. Since neutrino data requires at least one neutrino mass 
to be larger than about 0.05 eV, we do not find any solutions for My > A2 x 10 15 GeV and 
M24 > 8x 10 14 GeV, respectively. In the latter case we have assumed that all 24-plets have 
similar masses. For sake of completeness we note that one can also satisfy all neutrino data 
by giving one of the 24-plets a large mass in the order of Mqut or larger having a model 
with effectively only two 24-plets. 



A. Lepton flavour violation 

We have seen in eq. ( 1241) that rates for the lepton flavour violating decays of \x and r scale 
like the LFV entries in the slepton mass squared matrix squared and inverse to the overall 
SUSY mass to the power eight. From this one immediately concludes the rates for the rare 
lepton decays are in general larger in seesaw models of type-II and III than in type-I models 
for fixed SUSY masses and seesaw scales except if one arranges for special cancellations. 

Comparing the type-II with the type-Ill model one finds that LFV decays are larger for 
type-Ill, as shown for the case of fi — > e'j in Fig. [6j From eqs. (1211 and (l22l) . however, one 
would expect that type-II should have larger LFV. Numerically we find the opposite for two 
reasons, (i) Briji — > Ijj) strongly depends on the SUSY masses, see eq. (124)) and type-Ill 
has a lighter spectrum than type-II (for the same mSUGRA input parameters). And (ii) 
2-loop effects are very important in type-Ill, due to the large coefficients, in general leading 
to large flavor violating soft SUSY breaking parameters. 

In Fig. [6] we compare Br(li — > Ijj) for the three seesaw models taking degenerate seesaw 
spectra in case of type-I and type-Ill. Note that in case of seesaw type-Ill we can only 
show a relatively short interval for the seesaw scale which is mainly due to two reasons: 
(i) for scales below approximately 10 13 GeV the gauge couplings get large at Mqvt as a 
consequence of the large beta functions and, thus, perturbation theory breaks down, (ii) 
One encounters negative mass squares for the scalars, in particular for the lighter stau and/or 
lighter sbottom. The latter point is also the reason why the possible range is larger in case 
of the larger soft SUSY breaking parameters. 



The values for Br(fi — > cy) in Fig. [6] are larger than the current experimental bound [58 
so one might worry if in case of type-Ill models only SUSY spectra beyond the reach of 
the LHC are allowed. (Note, that even for the examples shown the masses of the sfermions 
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m =M 1/2 =300 (GeV), tan|3=10, A =0 (GeV) m =M 1/2 =1000 (GeV), tanp=10, A =0 (GeV) 




FIG. 6: Br(fj> — > cy) as a function of the seesaw scale for seesaw type-I (red line), seesaw type-II 
(blue line) and seesaw type-Ill (magenta line). In case of type-I and type-Ill a degenerate spectrum 
has been assumed. On the left panel mo = w-i/2 = 300 (GeV), on the right panel mo = w-i/2 = 1000 
(GeV). In both cases we take tan/3 = 10, Aq— and [i > 0. 

m =M 1/2 =1000 (GeV), tanp=10, A =0 (GeV) 




s 13 

FIG. 7: Br(n ->■ ej) versus sj 3 for m = M l/2 = 1000 GeV, tan /5 = 10, A = GeV and /x > 0, for 
seesaw type-I (solid lines) and seesaw type-Ill (dashed lines), for Ms ecsaw = 10 14 GeV. The curves 
shown are for 2 values of the Dirac phase: 5 = (red) and 5 = tt (blue), both for normal hierarchy. 

are already in the range of several hundred GeVs as can be seen from tabled) Indeed we 
find that by putting generic Yukawa couplings which are able to explain neutrino data one 
needs a heavy spectrum to be consistent with bounds on the rare lepton decays. However, 
this is strictly true only for the TBM angles and R =1. Accidental cancellations due to 
different contributions to the flavor violating soft masses and thus to the rare lepton decays 
are possible in type-Ill (and in type-I). As an example we show in Fig. [7]i?r(/i — > ey) as a 
function of the reactor angle sf 3 for different values of the Dirac phase 5. For comparison 
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TABLE I: Examples masses in GeV for M 1/2 = 1000 GeV, tan /3 = 10, A = GeV and /i > 0, for 
seesaw type-Ill for a degenerate seesaw spectrum with M24 = 10 14 GeV. 

m =M 1/2 =1000 (GeV), tanp=10, A =0 (GeV) 
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FIG. 8: Br{r — > erf) versus (left) and Br(r — > fi'y) versus sf 3 (right) for too = M X j 2 = 1000 
GeV, tan (3 = 10, Aq = GeV and /i > 0, for seesaw type-I (solid lines) and seesaw type-Ill (dashed 
lines), for Mg eesaw = 10 14 GeV. The curves shown are for 5 = (red) and 5 = tt (blue) for normal 
hierarchy. 

we also show the calculation for a type-I model. For 5 = n there is a range of s^ 3 where this 
branching ratio is below the experimental constraint. 

At first glance this seem to require some fine-tuning of the underlying parameters. How- 
ever, one can look at this from a different perspective: Assume that the MEG collaboration 
has found a non- vanishing value for Br(fi — > e^j) and from LHC data one has found that the 
spectrum is consistent with the type-Ill seesaw model. For a fixed i?-matrix, e.g. R=l one 
would obtain in this case a relation between and M24. This can be exploited to put a 
bound on M24 or even to determine it depending on the outcome of measurements of reactor 
angle and, thus, the model assumptions can be tested. In Fig. [8] we show the corresponding 
rare tau decays. Note that also for r — e7 such a cancellation exists in principle but the 
corresponding range is excluded by fi — > e^y. In contrast r — > /ry is insensitive to the reactor 
angle and should be measurable in the near future. 

Up to now we have assumed that the seesaw spectrum is nearly degenerate which is of 
course a strong assumption. We show in Fig. |9]two examples where we keep in each case two 
masses fixed and vary the third one. Note, that in contrast to SUSY particles the indices 
of the heavy particles are generation indices and do not correspond to a particular mass 
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FIG. 9: Branching ratios for — > Ijj (solid lines) and Zj — ?• (dashed lines) versus the seesaw scale 
for tan j3 = 10, \i > 0, = GeV, M1/2 = mo = 1000 GeV. On the left panel we scan on Mrj with 
M/j 2 = Mr 3 = 2 x 10 13 GeV while on the right panel we scan on Mr 3 with = M/j 2 = 2 x 10 
GeV. The color code is red for fi — >■ or fj, — > 3e, blue for r — )■ fMj or r — >■ 3// and green for r — ?■ e7 
or r — >• 3e. 

ordering, e.g. M/j 2 corresponds to the 'solar neutrino scale' and Mr 3 to the 'atmospheric 
neutrino scale'. In case that the mass of the first generation state is varied, e.g. the left plot 
of this figure, one finds a decrease of the branching ratios with increasing seesaw mass . 
This is mainly caused by an increase of the SUSY spectrum while at the same time neutrino 
physics is only affected mildly requiring only a light increase of the corresponding Yukawa 
couplings to obtain the correct neutrino masses. If, on the other hand, the mass Mr. 3 of the 
third generation seesaw particles is increased on needs also a sizable increase of the Yukawa 
couplings to obtain the correct neutrino mass difference squared for the atmospheric sector. 
This leads to the observed behaviour that the branching ratios for r — > ^7 and r — > 3fi 
increases while the other ones decrease. 



B. Dark Matter 

The changes in the spectrum induced by the new heavy states also impact on the pre- 
dictions with^espect to the relic density which we have calculated using the program 
micrOMEGAs 
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As is well-known, within mSUGRA there are 4 regions in parameter 
space, in which the constraint from dark matter can be satisfied. These are (i) the bulk re- 
gion; (ii) the stau co-annihilation region; (iii) the focus point line and (iv) the Higgs funnel. 



Below we will show usually the range of Qh 2 allowed at 3 a according to [58j 

0.081 < Qh 2 < 0.129. (25) 
In particular, the co-annihilation region is very sensitive to the difference between the 
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FIG. 10: Difference between the masses and the lightest stau and the lightest neutralino (upper 
row) as well as the corresponding flh 2 (lower row) as a function of the seesaw scale. The left (right) 
plots are for seesaw type-II (III). A degenerate seesaw spectrum has been assumed in case of seesaw 
type-Ill. Mi/2 = 800 GeV, Aq = 0, tan/3 = 10 and [x > 0. The lines correspond to full blue line 
mo = 0, red dashed line ttiq = 50 GeV, green dashed dotted line tuq = 100 GeV, black dashed 
line mo = 150 GeV and orange full line tuq = 200 GeV. The gray band shows the preferred range 
according to eq. (f25|) . 

masses of the lightest stau and the lightest neutralino. In Fig. [10] we observe that this differ- 
ence depends strongly on the seesaw scale in both models. For a fixed M1/2 and m lowering 
the seesaw scale increases this mass difference, which then leads to a larger calculated Qh 2 . 
To compensate for this effect one needs to lower m , with the value depending on the seesaw 
scale chosen. For certain seesaw scales then tuq needs to be lowered below tuq = and the 
co-annihilation region disappears. In this region of parameter space both models behave in 
a qualitatively similar way. However, recall that spectra run faster towards smaller masses 
in seesaw type-Ill. 

Also the focus point region is very sensitive to the precise values of the input parameters. 
The focus point region appears in mSUGRA for large values of m and small/moderate 
values of M1/2 of the order of (9(100) GeV, the exact value depending on m®. This can be 
seen in figs. [TT1 and [T2l where we show m^o, the higgsino content 1 7V"i 3 1 2 + | A r i4 1 2 and the 
corresponding Qh 2 as a function of tuq for a fixed seesaw scale Mt,w — 10 14 GeV, Aq = 0, 
tan/3 = 10, fi > and various values of Mi/2- Note, that we take different values of Mi/ 2 
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FIG. 11: Mass of the lightest neutralino (left plot), its higgsino content (middle plot) and the 
corresponding O/i 2 (right plot) as a function of ttiq for a seesaw type-II model with My = 10 14 
GeV, rritop = 171.2 GeV, Aq = 0, tan/3 = 10 and \x > 0. The lines correspond to full blue line 
M 1/2 = 195 GeV, red dashed line M 1/2 = 200 GeV, green dashed dotted line M 1/2 = 205 GeV, 
black dashed line M]y 2 = 210 GeV and orange full line Mi/ 2 = 215 GeV. The gray band shows the 
range eq. ([25]) . 
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FIG. 12: Mass of the lightest neutralino (left plot), its higgsino content (middle plot) and the 
corresponding fth 2 (right plot) as a function of mo for a seesaw type-Ill model with a degenerate 
seesaw scale My/ = 10 14 GeV, mt op = 171.2 GeV, Aq = 0, tan/3 = 10 and /i > 0. The lines 
correspond to full blue line M-^m = 400 GeV, red dashed line M\j 2 = 405 GeV, green dashed dotted 
line Mi/2 = 410 GeV, black dashed line M 1/2 = 415 GeV and orange full line M 1/2 = 420 GeV. 
The gray band shows the range eq. (125 p . 

for the two models in such a way that we obtain similar values for m^o. We find that 
both models behave differently in this region of parameter space, e.g. the higgsino content 
|-^i3| 2 + |-^i4| 2 decreases (increases) with increasing values itlq for seesaw type-II (type-Ill). 
However, also for type-II the higgsino content increases for increasing rriQ once we reach 
the multi-TeV range but we did not get correct electroweak symmetry breaking in case of 
multi-TeV values for m in case of type-Ill models. The increased higgsino content of the 
lightest neutralino leads to on increase (decrease) of its couplings to the Z-boson and the 
light Higgs boson (to sfermions) resulting in the observed dependence of Qh 2 for tjiq close 
to the 1-TeV region. 

With these observations it is clear that the DM allowed regions will be shifted in the 
rriQ-Mi/2 plane compared to the usual mSUGRA expectations. We fix in the following 
m top = 171.2 GeV, tan/3 = 10, A = and \x > as well as the seesaw scale to 10 14 GeV. 
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FIG. 13: Dark matter allowed region (in blue) for mSUGRA (left panel) and for type-I seesaw 
(right panel). The parameters are tan/3 = 10, Aq = 0, [i > and Mt = 10 14 GeV for mtop = 171.2 
GeV. Also shown (in yellow) are the regions excluded by LEP (small values of M1/2), and by LSP 
constraint (small values of mo). Also shown are the Higgs boson mass curves for Mh = 110 GeV 
(in red) and for = 114.4 GeV (in magenta). 



For comparison we show in Fig.[T3]the usual mSUGRA case without any heavy intermediate 
particles (left plot) as well as the case of a seesaw type-I scenario (right plot). The blue 
bands show the 3er range according to |58| and we see the three usual regions: the stau 
co-annihilation with a lighter stau mass close to the LSP mass for M1/2 < 300 GeV, the 
bulk region for moderate values of M1/2 and m resulting in small sfermion masses as well as 
the focus point region for Mi/ 2 — 170 GeV and large values of m . In addition, we show the 
lines corresponding to Mh = 110 GeV and 114 GeV. Note, that the theoretical uncertainty 
on Mh is still of the order of 3-5 GeV 
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l61| . Moreover, the value of the Higgs boson 
mass also depends strongly on Aq and in particular for negative values of Aq one can easily 
increase the value of Mh while the DM allowed regions hardly change. 

The part of parameter space most affected is the one at large m . Since in mSUGRA 
fj, is calculated from the requirement of correct electroweak symmetry breaking, \x changes 
rapidly in this region. With the Higgsino content in the lightest neutralino changing rapidly 
as a function of /i, this region is then very sensitive to any changes of parameters. Since 
the Y v also impacts on the running of the Higgs mass parameters and thus slightly affects 
the value predicted for /x, some small changes are found relative to mSUGRA here. Note, 
however, that this region is hi ghly constrained by the lower bound on the lightest chargino 
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mass of the order of 103 GeV 

In case of the other two seesaw models the shift of the allowed regions is much more 
pronounced, as discussed above. In Figs. [T3] and [15] we show to regions for type-II (left 
plot) and type-Ill (right plot) and two different values for Aq. As claimed above, the Higgs 
mass bounds gets shifted significantly while the DM allowed regions are hardly affected. As 
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FIG. 14: Like in Fig. [13] but for seesaw type-II (left panel) and type-Ill (right panel). 



M T = 10 (GeV) tanp=10, A =-300 (GeV) 
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FIG. 15: Like in Fig. [TJ]but for Aq = —300. Seesaw type-II (left panel) and type-Ill (right panel). 



expected the effects are much more pronounced in case of type-Ill as the effects of the heavy 
particles on the spectrum is much stronger. Note, that in particular the bending of the 
allowed region for large m is due to the changed higgsino content as discussed in case of 
figs. [11] and [12] Moreover, the case of stau co-annihilation is not viable anymore in case of 
the type-Ill model already for this value of the seesaw scale. For completeness we mention 
that for the type-II the stau co-annihilation region disappears (below Mu^ = 1500 GeV) for 
Mt %j 10 13 GeV. For completeness we note that the results here differ slightly from the ones 
of our previous work 42| because (i) of the corrections of the 1-loop RGEs of ref. 14| by 



47j and (ii) the complete set of 2-loop RGEs are now used. 



In the case of large tan f3 an additional region, usually called the Higgs funnel, opens up. 
This region is characterized by ~ 2m^o. Also here the regions gets shifted compared to 
usual mSUGRA scenario. However, this region is very sensitive to higher order corrections 
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FIG. 16: Comparison between using 1-loop (red) or 2-loop (blue) RGEs on the dark matter allowed 
region for type-II (left panel) and type-Ill (right panel). The parameters are: Aq = 0, fi > and 
M S ee S aw = 10 14 GeV, m top = 171.2 GeV and tan/? = 52 for type-II and tan/3 = 49 for type-Ill. 



and therefore it is quite important to use full 2-loop RGEs as can be seen in Fig. [T6j We 
have again fixed Aq = 0, /x > 0, m t0 p = 171.2 GeV and the seesaw scale to 10 14 GeV, with 
a degenerate spectrum in case of the type-Ill model. The main reason for the observed 
and rather surprisingly large differences between the different calculations is that the 2-loop 
contributions decrease the neutralino mass compared to the 1-loop case while at the same 
time increasing M4. For example, in case of seesaw II and for fixed values of tuq = Mi/ 2 = 
1500 GeV we get in case of 1-loop RGEs m^o = 560 GeV, M4 = 1090 GeV and in case of 
2-loop RGEs m^o = 498 GeV, Ma = 1100 GeV. For completeness we note that this region 
is also very sensitive to input values for m t and mj, 42 



IV. CONCLUSIONS AND OUTLOOK 



To summarize, we have investigated in detail a supersymmetric version of a seesaw model 
of type-Ill and compared it to seesaw models of type-I and type-II. In case of type-II and 
type-Ill models we have embedded the £77(2) triplets in the corresponding SU (5) represen- 
tations to maintain gauge coupling unification, e.g. 15-plets in case of type-II and 24-plets 
in case of type-Ill models. For definiteness we have assumed mSUGRA boundary conditions 
for the soft SUSY breaking parameters. 

The additional heavy charged states lead to changes in the beta-functions and, thus, also 
in the running of the SUSY mass parameters. We have calculated the soft masses as a func- 
tion of the seesaw parameters. As discussed in some detail, there are certain combinations 
of soft masses, which are approximately constants over large regions of mSUGRA space. 
These "invariants" contain indirect information about the seesaw scale assuming the type of 
seesaw model. In certain parts of the parameter space, e.g. for low seesaw scales, one might 
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even be able to exclude certain seesaw models by combining mass measurements at the LHC 
with the mSUGRA paradigm. We note, that using 2-loop RGEs will be crucial to obtain 
reliable results. 

The changes in the spectrum leads obviously to changes in the phenomenology. We have 
calculated lepton flavour violating observables, such as Br(l{ — > lj + 7). We find that for 
fixed (degenerate) seesaw scale these branching ratios are in general largest for type-Ill 
models followed by type-II and type-I. This is a consequence of the fact that for a given set 
of mSUGRA parameters the spectrum in type-Ill is lighter than for type-II models which 
is again lighter than in type-I models. However, the difference in the predictions of type-II 
and type-Ill is somewhat smaller than expected from these considerations because in type-II 
models the flavour violating entries are larger compared to the case of type-Ill models. 

We also investigated the predictions for the relic density Qh 2 in the type-Ill model and 
compared them with the other models. We find the usual four regions in the mSUGRA 
parameter space but of course they are shifted due to the changes in the spectrum. It has 
been found that in particular in case of the Higgs-funnel the use of 2-loop RGEs is crucial to 
identify the correct allowed region. Last but not least we note, that for low seesaw scales the 
co-annihilation region vanishes for both, the type-II and the type-Ill models, as the required 
mass difference between the lightest neutralino and the stau cannot be obtained anymore. 
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Appendix A: RGEs for the seesaw type-II and seesaw type-Ill models at 2-loop 

In the appendix we collect the beta coefficients for the gauge couplings as well as anoma- 
lous dimensions of the superfields which are the ingredients to calculate the 2-loop RGEs 
for both, the seesaw type-II and type-Ill, models using the procedure given in |63| based 
on the spurion formahsm & The complete set of RGEs for both models at 2-lo^ is also 



given online 69J. In the following we briefly summarize the basic ideas of this calculation 



for completeness. 
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For a general N = 1 supersymmetric gauge theory with superpotential 
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The /3-functions for the superpotential parameters can be obtained by using superfield 
technique 65|, |66j. The obtained results are 46| 
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The exact results for the soft /3-functions are given by [63] : 
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where 7 is the matter multiplet anomalous dimension, /3 g the beta function for the gauge 
coupling g; the (..) in the superscripts denote symmetrisation and 
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Here M is the gaugino mass and Y ZJ 
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(1A6I) hold in a class of renormalisation schemes that includes the DRED'-one 67fl. We take 
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C(R), C(G) being the quadratic Casimirs for the matter and adjoint representations, respec- 
tively, Q = T(R) — 3C(G), and rT(R) = tr [C(i?)], r being the number of group generators. 

In the following subsections we give the anomalous dimensions and beta functions for the 
seesaw models of type-II and type-Ill so that with the help of the above equations one can 
calculate all RGEs at the 2-loop level. 
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1. The anomalous dimensions for seesaw type-II 



Here and in the subsequent sections 1 denotes the 3x3 unit matrix. Nx is the number 
of generations of heavy field X. Furthermore, we define 
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= J^9 3 i{lQN s + 18N T + 66 + N z 
1 



(A41) 



150 



1194s 2 + 810^2 + 2640s 2 - 810(| Ai| 2 + |A 2 | 2 ) + (256s 2 + 3200s 2 ) N s 



(648s? + 2160S2 2 ) N T + (s? + 45s 2 2 + 80s 3 2 ) N z - A20Ti[Y d Y c 
540Trfr e F e tN ) - 720Tr(V s y;) - 810Tr(Y t y t *) - 780Tr (y u Y^\ - 420Tr(yyj' 



1 



^l ) =^gl[SN z + 4N T + 2 

P$ = ^2 [l8Si 2 + 250s 2 2 + 240s 3 2 - 70|A!| 2 - 70|A 2 | 2 + (48s? + 240s 2 2 ) N T 

+ (s 2 + 105s 2 2 + 80s, 2 ) N z - 60IY {Y d Y^ - 20Tr(Y e F e t) - 70Tr(yy* 

- 6(m(y u y u t) - 60Tr(Y z y;) 



PS 
0£ 



3 + N z )+5N l 



15 



33s? + 135s 2 2 + 210s 3 2 + 5 ( 145s 2 + 8^ ) N s - 135Tr [Y^ 



(A42) 
(A43) 



(A44) 
(A45) 



26 



+ (g{ + Ahg 2 2 + 170^|) N z - 60Tr (yyj) - 601^1^) - 60Tr (y z Y^ 



(A46) 



The anomalous dimensions for seesaw type-Ill 



^ (45g 2 2 + SOgj + gf) 1 + YjY d + Y^Y U (A47) 



30 

\g\Y^Y u - TY\Y d Y\Y A - 2Y^YjY*Y d - 2Y^Y U + 1 



199 4 

900^ + T0 9 ^ + ~4 92 



+ ^ghl + - ^3 + % 3 4 AW + 3g 4 2 N WM + [1^ + ^ 2 4 + |g ) iv Xjv/ 



- ^y u Tr(y 6 yj) + yjy ( - 3Tr(yyj) + - Tr(yy e t)) - 3Y^Y u Tt(yX 

- hjY u Tr (y^) - 3YjY u Tr (y.Yj) (A48) 
^ (l0y e ty + lSY^ - 3 (bg 2 2 + <7 2 ) 1 + 31?n) (A49) 

[24o^Yjy e + 1200^^ - sey^yyjy - eoy/yyjy, - 4ooy e t yy e t y e 

- 45Yj;yy b t y - 300YiyyJy w + eif ioo<? 2 4 av m + 125^ + 25(3^ + # 4 ) a> Xm 



+ 3o^ 2 2 + 69s 4 ) - isy^n Tr (yyj) - goy^y, Tr (yy/) - eooyjy Tr (yyj) 

- 2ooy e t y e Tr(y e y e t) - isoy/yTr (yyj) -9ooyjy„Tr(yyj) -9oy 6 f yTr(yyJ 

- 45oyiy w Tr (y w y^) -isoy/yTr(y y;) -9ooy^Tr(yyi) ] (aso) 

3 Tr (y d yj) - 1 (s^ 2 + <? 2 ) + Tr (Y e Y^) ( A51) 



- ^ 2 Tr(yy;) + i6s 2 Tr(yy]) + ^ 2 Tr(yy e t) _ 1 Tr (yy^yy/) 

- 6Tr(yy]yJy;) - 3Tr(y d y u ty 1 yj) - 3 Tr (y e y e ty e y e t) -^Tr (yy^yj 



(A52) 

10 Tr (y u y u t) - 10 Tr (y x Y^) + %2 " 5 Tr (y^) - Tr (y Yj) + g{ 



10 



+ + ^ 9 bl + ^at + 3^ + I (3s 2 4 + £ 4 ) N Xm + \g\ Tr (y 

+ 16^ 2 Tr (yyj) + 6s 2 Ti (Y w Y^j +2g 2 Ti(Y x Y^ + 16^ 2 Tr (y x Y}) 

- A ^ (yy e f y n f ) - ^ Tr (yyiyy/) - e Tr (y y] yJy;) - 3 Tr (y y! yyi) 



(A53) 



27 



- 1 Tr {Y e YlY w Y}) - 9 Tr {y u Y^Y u Y^) - H Tr (y^y^) - 9 Tr (V^Fj) 
27 / \ 

-— Tv(y 6 y 6 ty 6 y 6 tj (A54) 

2 n + >7^ T ) " ^ (20^| + <??) 1 (A55) 



+ IglY^Yj + 6g 2 2 Y d *Y d T - 2Y^Y X - 2Y d *Y d T Y d *Y d T - TY^Y^J 



225 



imglgl + 1800s 3 4 7V Gm - 2(% 3 4 + 202,? 4 + 75 (8g 4 3 + gf)N 



*x M 



- QY d *Yj Tr (y d yj) - 2Y*Yj Tr (y e y e t) + Y*Y X (lg\ - 3 Tr (y w Yi) + 6gl 

- 6 Tr (y u yj) - 6 Tr (y x Y^) " Tr faf 



(A56) 
(A57) 



- 2{y:yJy:yJ + y;y u T y;y„ T ) + — 1 [i6o<& 3 2 + 2i% 4 + 450s 3 4 a Gm - 5o^ 4 

1. 



75 [2gt + g{ )N 



— Y*Y T 
5 u " 



15 Tr (Y W YS) + 2# 2 - 30s 2 + 30 Tr {y u Y^ 



+ 30 Tr (Y X YA +3Tr( Y b Y 



2Y e *Y e T - Q -gll 

^ V*vT\r*\rT 
^ x e 1 b I b I e 



2y;y c t y;y c t - 3y;yJX*y e T + ^1 (2bn Xm + 78 



+ Y*Y T 

' e e 



2 Tr ( y e y e t) + Qg 2 2 - 6 Tr (Y d Yj 



> 2 

-9i 



-% 2 i + y w x T 



2g A 2 l(3N XM +AN WA 



(A58) 
(A59) 

(A60) 
(A61) 



101 + Io 



w e e w 



1 x\s*Y T Y*Y T 

" t» IB W W 



yjyjj - ios 2 - i5Tt(y w y^ - 30Tr(y u yj) _ 3 oTr(yyj- 



3Tr (Y b YA +6gl 



6gtl(3N GM 

3, 



3 + N 



rY b *Y b 



T 



(A62) 
(A63) 
(A64) 

(A65) 



3 

50 



oy*yTy*yT r / 9 y* yTy* yT , q v* V r V* V 
6I b *b I b 1 b ~ °l AI b 1 e 1 e 1 b "> 0I b *w I w 1 l 

3y;y fe T {ios 2 -ioTr(yy u t) -ioTr(yyj) +2 5 2 -5Tr(y^) - Tr(y 6 y f )} 

(A66) 

i(l6^ 2 + 5^ 2 + 9^)l (A67) 



28 



(2) 1 -, 

7 - = — 1 

'x M 3 6 



'x M 6 



223^ + 90g 2 l9 2 2 + 135^| + 160g 2 l9 2 3 + 288(g 2 2 g 2 3 + «$V G J - 32^ + 108<7 4 AV M 
(75^ + 81gl + 96^|) A> Xm ] (A68) 
l - ( - (165 2 + 5gf + % 2 2 ) 1 + 6F;Fj) (A69) 



7 



(2) 



+ — 1 

36 



223^ + Wgbl + 135^| + \mg\gl + 288(^ 3 2 + <#V G J - 32^ + lWg A 2 N WM 

20 fc^yjyj + Y x *Y^Y x *Yf 



10 



+ (75(7 4 + 81^| + 96s 4 ) N Xl 
-Y^{l5T*(Y„Yl)+30Tr(^ + Ag\ 



(A70) 



4. Beta coefficients for the seesaw type-Ill at 2-loop level 



1 



66 



(A71) 



150 



6(199^ + 135#| + 440^1 - 9 Tr \Y b Y^ 
70 Tr (Y d Yf) - 90 Tr (V e Y e +) - 130 Tr (y u Y?) - 60 Tr (y w Y^) - 190 Tr (V^) } 



^=^2(3^+4^ + 2 



(A72) 
(A73) 



30^ 



54(?f + 750s 2 + 720(4 2 + g 2 2 N WM ) + 15 I64 2 + 21^ + 5gf N x 



l8Tt(Y b Y^ - 280Tr (y w Y^ - 60 Tr (Y e Y e ^ - 180 (Tr (Y d Y^j + Tr (y u Y^ 
Tr(y x F x t))] (A74) 

(A75) 



33^ + 135^ + 210^ + 810^A^ Gm + 5 34^ + 5g( + 9gZ)N XM 



170giN XM -60Tr(Y d Yj 



60Tr(y u yJ 



eoTrfyyJ 



(A76) 



[1] Y. Fukuda et al. [Super-Kamiokande Collaboration], Phys. Rev. Lett. 81, 1562 (1998) 
[2] SNO, Q. R. Ahmad et al, Phys. Rev. Lett. 89, 011301 (2002), [nucl-ex/0204008] . 
[3] KamLAND, K. Eguchi et al, Phys. Rev. Lett. 90, 021802 (2003), [hep-ex/0212"021 . 



29 



[4] KamLAND Collaboration. larXiv:0801.4589l [hep-ex], 

[5] T. Schwetz, M. A. Tortola and J. W. F. Valle, New J. Phys. 10, 113011 (2008) [arXiv: 0808 .20161 

[hep-ph]]. Version 3 on the arXive is updated with data until Feb 2010 
[6] S. Weinberg, Phys. Rev. Lett. 43, 1566 (1979); S. Weinberg, Phys. Rev. D 22, 1694 (1980). 
[7] E. Ma, Phys. Rev. Lett. 81, 1171 (1998) [arXiv:hep-ph/9805219| . 
[8] P. Minkowski, Phys. Lett. B 67 (1977) 421. 

[9] T. Yanagida, in KEK lectures, ed. O. Sawada and A. Sugamoto, KEK, 1979; M Gell-Mann, 
P Ramond, R. Slansky, in Supergravity, ed. P. van Niewenhuizen and D. Freedman (North 
Holland, 1979); 

[10] R.N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44 912 (1980). 

[11] J. Schechter and J. W. F. Valle, Phys. Rev. D 22, 2227 (1980). 

[12] T. P. Cheng and L. F. Li, Phys. Rev. D 22, 2860 (1980). 

[13] R. Foot, H. Lew, X. G. He and G. C. Joshi, Z. Phys. C 44, 441 (1989). 

[14] A. Rossi, Phys. Rev. D 66, 075003 (2002) [arXiv:hep-ph/0207006] . 

[15] M. R. Buckley and H. Murayama, Phys. Rev. Lett. 97, 231801 (2006) [arXiv:hep-ph /0606088| . 
[16] J. Hisano, T. Moroi, K. Tobe and M. Yamaguchi, Phys. Rev. D53, 2442 (1996), 
[hep-ph/9510309| . 

[17] J. R. Ellis, J. Hisano, M. Raidal and Y. Shimizu, Phys. Rev. D 66, 115013 (2002) 
[arXiv:hep-ph/02061l"0] . 

[18] F. Deppisch, H. Paes, A. Redelbach, R. Riickl and Y. Shimizu, Eur. Phys. J. C28, 365 (2003), 
[hep-ph/0206122| . 

[19] S. T. Petcov, S. Profumo, Y. Takanishi and C. E. Yaguna, Nucl. Phys. B 676 (2004) 453 

[arXiv:hep-ph/0306195] . 
[20] E. Arganda and M. J. Herrero, Phys. Rev. D73, 055003 (2006), [| hep-ph/0510405] . 
[21] S. T. Petcov, T. Shindou and Y. Takanishi, Nucl. Phys. B 738, 219 (2006) 

[arXiv:hep-ph/0508243] . 

[22] S. Antusch, E. Arganda, M. J. Herrero and A. M. Teixeira, JHEP 11, 090 (2006), 
[hep-ph/0607263| . 

[23] F. Deppisch and J. W. F. Valle, Phys. Rev. D 72, 036001 (2005) [arXiv:hep-ph/0406040] . 
[24] J. Hisano, T. Moroi, K. Tobe, M. Yamaguchi and T. Yanagida, Phys. Lett. B357, 579 (1995), 
[hep-ph/9501407| . 

[25] M. Hirsch, J. W. F. Valle, W. Porod, J. C. Romao and A. Villanova del Moral, Phys. Rev. D 

78, 013006 (2008) [arXiv: 0804 .40721 [hep-ph]]. 
[26] E. Arganda, M. J. Herrero and A. M. Teixeira, JHEP 10, 104 (2007), [0707.2955]. 
[27] F. Deppisch, T. S. Kosmas and J. W. F. Valle, Nucl. Phys. B752, 80 (2006), [hep-ph/05 12360] , 
[28] M. Hirsch, S. Kaneko and W. Porod, Phys. Rev. D 78, 093004 (2008) |arXiv:0806.336T1 [hep- 

Ph]]- 

[29] C. Biggio and L. Calibbi. ia7Xiv:1007.3750l [hep-ph], 

[30] J. Hisano, M. M. Nojiri, Y. Shimizu and M. Tanaka, Phys. Rev. D 60, 055008 (1999) 
[arXiv:hep-ph/98084T0] . 



30 



[31] J. N. Esteves, J. C. Romao, A. Villanova del Moral, M. Hirsch, J. W. F. Valle and W. Porod, 

JHEP 0905, 003 (2009) [arXiv: 0903 .14081 [hep-ph]]. 
[32] G. A. Blair, W. Porod and P. M. Zerwas, Eur. Phys. J. C27, 263 (2003), |hep-ph/0210058 . 
[33] A. Freitas, W. Porod and P. M. Zerwas, Phys. Rev. D72, 115002 (2005), [hep-ph/0509056] . 
[34] F. Deppisch, A. Freitas, W. Porod and P. M. Zerwas, 0712.0361. 

[35] B. C. Allanach, J. P. Conlon and C. G. Lester, Phys. Rev. D 77, 076006 (2008) jarXiv: 080 1.3666 1 
[hep-ph]]. 

[36] A. Abada, A. J. R. Figueiredo, J. C. Romao and A. M. Teixeira. larXiv: 1007.48331 [hep-ph], 

[37] E. Komatsu et al.. larXiv: 1001. 45381 [astro-ph.CO]. 

[38] K. Griest and D. Seckel, Phys. Rev. D 43, 3191 (1991). 

[39] H. Baer, C. Balazs and A. Belyaev, JHEP 0203, 042 (2002) [arXiv:hep-ph/0202076| . 

[40] J. L. Feng, K. T. Matchev and T. Moroi, Phys. Rev. D 61, 075005 (2000) 

[arXiv:hep-ph/9909334] . 
[41] M. Drees and M. M. Nojiri, Phys. Rev. D 47, 376 (1993) [arXiv:hep-ph/9207234| . 
[42] J. N. Esteves, S. Kaneko, J. C. Romao, M. Hirsch and W. Porod, Phys. Rev. D 80, 095003 

(2009) |arXiv:0907.5090l [hep-ph]]. 
[43] F. Staub. Ia7xiv:0806.0538l [hep-ph], 

[44] F. Staub, Comput. Phys. Commun. 181, 1077 (2010) [arXiv: 0909. 28631 [hep-ph]]. 
[45] F. Staub. laTXiv: 1002.08401 [Unknown], 

[46] S. P. Martin and M. T. Vaughn, Phys. Rev. D 50, 2282 (1994) [Erratum-ibid. D 78, 039903 

(2008)] | arXiv:hep-ph/9311340 . 
[47] F. Borzumati and T. Yamashita. larXiv: 0903. 27931 [hep-ph], 
[48] S. K. Kang, T. Morozumi and N. Yokozaki. larXiv:1005.13"5l [hep-ph], 
[49] J. A. Casas and A. Ibarra, Nucl. Phys. B 618, 171 (2001) [arXiv:hep-ph/0103065| . 
[50] L. Calibbi, M. Frigerio, S. Lavignac and A. Romanino, JHEP 0912 (2009) 057 |arXiv:0910.0377l 

[hep-ph]]. 

[51] L. Calibbi, Y. Mambrini, S. K. Vempati, JHEP 0709 (2007) 081. |arXiv:0704.35l8l [hep-ph]]. 
[52] B. C. Allanach et al., in Proc. of the APS/DPF/DPB Summer Study on the Future of Particle 

Physics (Snowmass 2001) ed. N. Graf, Eur. Phys. J. C 25 (2002) 113 [arXiv:hep-ph/0202233 . 
[53] W. Martens, L. Mihaila, J. Salomon and M. Steinhauser, arXiv: 1008.3070 [hep-ph]. 
[54] W. Porod, Comput. Phys. Commun. 153, 275 (2003) [arXiv:hep-ph/0301101| . 
[55] For the latetst version of SPheno, see the web page: 

http:/ /www.physik. uni-wuerzburg.de/~porod/SPheno. html 
[56] S. Antusch and M. Ratz, JHEP 0207, 059 (2002) [arXiv:hep-ph/0203027| . 
[57] P. F. Harrison, D. H. Perkins and W. G. Scott, Phys. Lett. B 530, 167 (2002) 

[arXiv:hep-ph/0202074] . 
[58] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667, 1 (2008). 

[59] G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, Comput. Phys. Commun. 176, 367 

(2007) | arXiv:hep-ph/0607059| . 
[60] B. C. Allanach, A. Djouadi, J. L. Kneur, W. Porod and P. Slavich, JHEP 0409, 044 (2004) 



31 



|arXiv:hep-ph/0406166] . 

[61] S. Heinemeyer, W. Hollik, H. Rzehak and G. Weiglein, Eur. Phys. J. C 39, 465 (2005) 



|arXiv:hep-ph/0411114 



[62] LEPSUSYWG, ALEPH, DELPHI, L3 and OPAL experiments, note LEPSUSYWG/01-03.1 
( [http:/ /lepsusy.web.cern.ch/lepsusy/Welcome.html ). 



[63] I. Jack, D. R. T. Jones and A. Pickering, Phys. Lett. B 426, 73 (1998) [arXiv:hep-ph/9712542 
[64] Y. Yamada, Phys. Rev. D 50 (1994) 3537 [arXiv:hep-ph/9401241| . 
[65] P. C. West, Phys. Lett. B 137 (1984) 371. 

[66] D. R. T. Jones and L. Mezincescu, Phys. Lett. B 138, 293 (1984). 

[67] I. Jack, D. R. T. Jones, S. P. Martin, M. T. Vaughn and Y. Yamada, Phys. Rev. D 50 (1994) 
5481 [arXiv:hep-ph/9407291]. 



[68] I. Jack, D. R. T. Jones and A. Pickering, Phys. Lett. B 432 (1998) 114 [arXiv:hep-ph/9803405 



[69] http:/ /theorie. physik.uni-wuerzburg.de/~fnstaub/supplementary.html 



32 



